Generation of entangled photons in graphene in a strong magnetic field 
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Entangled photon states attract tremendous interest as the most vivid manifestation of nonlocality 
of quantum mechanics and also for emerging applications in quantum information. Here we propose 
a mechanism of generation of polarization-entangled photons, which is based on the nonlinear optical 
interaction (four- wave mixing) in graphene placed in a magnetic field. Unique properties of quantized 
electron states in a magnetized graphene and optical selection rules near the Dirac point give rise 
to a giant optical nonlinearity and a high rate of photon production in the mid/far-infrared range. 

PACS numbers: 78.67.Wj, 42.65.Lm, 42.50.Ct 



To date, the most widely used method of generating 
entangled photons is based on the spontaneous paramet- 
ric down-conversion in a nonlinear crystal possessing a 
second-order nonlinearity [HE]. In this process, a pho- 
ton from a strong pump field at frequency uj p splits into 
two signal photons, uj p — ui\ +U2 which can be entangled 
in polarization, frequency, and wave vector. Entangle- 
ment in the polarization degree of freedom is the most 
convenient one for applications. Another way to generate 
quantum-correlated photons through a parametric non- 
linear optical process is spontaneous four-wave mixing in 
the optical fibers, in which two pump photons are con- 
verted into two signal photons, 2uj p — u)\ +ll>2, utilizing a 
third-order nonlinearity of silica [3] . This process is obvi- 
ously compatible with fiber communication technologies, 
although it does not directly lead to polarization entan- 
glement. In order to achieve the latter, one needs to use 
two pumps with crossed polarizations and apply addi- 
tional signal processing. In both nonlinear processes the 
photon pair production efficiency is very low. An alterna- 
tive approach utilizing the radiative decay of biexcitons 
in semiconductor quantum dots allows photon pairs 
to be generated on demand but requires cooling down to 
liquid helium temperatures. 

Graphene has unusual electronic and optical proper- 
ties stemming from linear, massless dispersion of elec- 
trons near the Dirac point and the chiral character 
of electron states [S]. Magnetooptical properties of 
graphene and thin graphite layers are particularly pecu- 
liar, showing multiple absorption peaks and unique se- 
lection rules for transitions between Landau levels jS}- 
I12j . Recent progress in growing high-quality epitaxial 
graphene and graphite with high room-temperature mo- 
bility and strong magnetooptical response attracted a lot 
of interest and paved the way to new applications in the 
infrared optics and photonics [T51 - fT5] . The time is ripe 
to explore the nonlinear and quantum optical properties 
of a magnetized graphene and their applications. We 
have recently shown that graphene placed in a magnetic 
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field possesses perhaps the highest infrared optical non- 
linearity among known materials llj. Here we argue that 
an extremely strong nonlinearity of graphene in combi- 
nation with its peculiar properties of the Landau lev- 
els open new avenues for generation of the nonclassical 
light states, in particular polarization-entangled photons. 
A similar mechanism of photon entanglement may exist 
in topological insulators where the surface states have a 
Dirac-cone dispersion and demonstrate similar properties 
of magneto-optical absorption. 




FIG. 1: Geometry of the proposed experiment. Two pump 
fields at frequencies ujhf and ujlf normally incident on a 
sheet of graphene placed in a magnetic field B generate entan- 
gled photons with opposite sense of the circular polarization. 

The proposed scheme is shown in Figs. 1 and 2. Here 
the energies of the Landau levels for electrons near the 
Dirac point are given by e n = sgn(n)hu c yj\n\ , where 
n = 0,±1,±2..., oj c = y/2v F /l c , v F « 10 8 cm/s the 
electron Fermi velocity, and l c — y^hc/eB the magnetic 
length. We assume that the graphene is biased or doped 
so that the Fermi level is between the states with n = 
-2 and n = -1, i.e. the state n = -2 is occupied and the 
states above are empty in the absence of pumping. Two 
incident strong pump fields at frequencies ojhf and u>lf 
resonant to the transitions from n = -2 to n = 1 and from 
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FIG. 2: Energy levels and optical transitions involved in reso- 
nant parametric generation of entangled photons in graphene. 
Left: Landau levels near the Dirac point superimposed on the 
linear electron dispersion without the magnetic field. Right: 
A scheme of the entangled photon generation process in the 
four-level system of LLs with energy quantum numbers n = 
-2, -1, 0, 1 that were renamed as states 1,2,3, and 4 for conve- 
nience of notation. Right-hand side (RHS) or left-hand side 
(LHS) circular polarization of light indicated on the figure 
corresponds to the allowed transitions. 



n = -2 to n = -1 respectively, generate two signal fields 
with opposite sense of the circular polarization at fre- 
quencies and that are close to resonance with 
transitions from n = -1 to and from n = to 1. Note 
that these transitions have the same energy. Therefore, 
the presence of the unshifted n = Landau level en- 
ables convenient entanglement in the polarization degree 
of freedom for two photons with nearly equal energies. 
All transition frequencies are easily tunable with a mag- 
netic field. 

The polarizations for the allowed transitions are indi- 
cated in Fig. 2. Here LHS and RHS denote left-hand 
and right-hand circularly polarized light with polariza- 
tion vectors in the (x,y) plane of the graphene defined as 
e( T ) = (x T iyo) respectively. Peculiar selection 
rules for graphene, A|n| = ±1 as opposed to An = ±1 
for electrons with usual parabolic dispersion, allow tran- 
sitions with a large change in the principal quantum 
number n, such as the transition from n = -2 to 1. 
The dipole matrix elements of the allowed transitions 
dmn ~ hevp/(e n — £m) grow fast (~ A) with increasing 
wavelength, and reach a large magnitude in the mid/far- 
infrared range; e.g. \d\/e — 13 nm for the n = to n = 1 
transition in the field of 1 Tesla (at wavelength of 34 /im). 
This enables an extremely high resonant third-order non- 
linearity [IT]. Note also that the states n = -1, 0, and 
1 have low population when the intensities of the opti- 
cal pumps are below saturation and Hlo c 3> fcgT. These 
factors lead to a high rate of photon generation and high 
signal to noise ratio of entangled photons. 

In order to determine the optimal conditions for entan- 
glement and the photon generation rate, we solve cou- 



pled equations for Heisenberg operators of the electron 
and signal photon fields, assuming that the strong pump 
fields are classical. Consider quasiparticles ("electrons") 
on Landau levels described by stationary states |m) and 
energy levels e m . After introducing creation and annihi- 
lation operators of an electron, aj„|0) = \m),a n \n) = |0), 
one can define a coordinate-dependent density operator, 

P mn {r,t) = ^^a]. n (t)a r:7n (t) (1) 
r i 

where the index j numerates individual electrons and the 
summation is carried over all electrons within a small vol- 
ume AV r in the vicinity of a point with radius-vector r. 
The density operator in Eq. ([I]) is normalized to the total 
electron density N(r) according to (9 E \ J2 m p mm \^ e) = 
iV(r) where \^>e) = ^Z m C m \m) is the wave function of 
the electron subsystem satisfying ^ m \C m \ 2 — 1. As- 
suming that the operators in different points of space 
commute with each other, the commutation relations be- 
come 

[p qp (r) , p mn (r')} = S(r-r')(p mp (r)6 qn -6 m pPq n (r)). (2) 

Using the density matrix defined in Eq. ([I]) , one can write 
the Heisenberg operator of any physical quantity x(r, t) 
as x = XnmPmnir^). In particular, the optical polar- 
ization is given by P(r,t) = d nm p mn where d mn is the 
dipole matrix element. 

The Heisenberg-Langevin equation for the density op- 
erator takes the form 



Pmn — ^_ \ h"mvPvn Pmvhvnj ~r~-^mji(Pmri) "r^mm (3) 

independently on whether a, operators obey the com- 
mutation relations for fermions or bosons. In Eq. ^ 
hnm = £n5 n m — d nm E(r, t) is the matrix element of the 
Hamiltonian operator H — h nm d^d m describing interac- 
tion with the electric field E(r,i) in the dipole approx- 
imation and R m n the relaxation operator, for which we 
will choose the simplest form R m ^n — ~lmnPmn- 

The Langevin noise operator F mn satisfies F mn = F^ m 
and (F mn ) = 0. Here the averaging (...) is taken both 
over the reservoir and over the initial state \^e) of the 
electron system. Its commutator is equal to 

[F mn (r,t),Fl n (T',t')] = 2 lmn (p nn -p mm )5{t-t')5(v-r') 

and correlation relations for its spectral components are 

(Ft mn (r')F u , mn (r)) = ^{p mm )5{u - J)5{v - r'), 



(^cj';mn( r )^ c I;mn( r )) 
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Their derivation follows standard steps described in 
|16) for any stationary random delta-correlated process, 
but without making the assumption of a thermal equi- 
librium. In the limiting case of a thermal equilib- 
rium, (p nn )/{p mm ) = exp (ftw nm /k B T), one can derive 
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from the above correlation relations that the fluctua- 
tion component of the polarization generated by the 
Langevin noise satisfies the fluctuation-dissipation the- 
orem [TBI [171. 

For a monochromatic electric field of a given field mode 
propagating in a dispersive medium with refraction index 

n(w), E = E e-' lut+lkz +~E,le~ luJt+lkz , one can introduce 
the operators of annihilation and creation of "photons 
in a medium" cq and c [TH] as E = eo-Eoc"o(r, t), Eg = 
eQ-EocJ(r, t). Here e is a unit vector of the polariza- 
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velocity. With this normalization of the field opera- 
tors the energy of the field in a volume V is given by 

W = Huj (^VcqCq + |J and their commutation relation 

reads [co,Cq] = y. Therefore, these operators determine 
the number density of the photons in a certain state \^f) 
of the field as nph = ('f f|cJco| , I'_f). This normalization 
is more convenient for field propagation problems than 
the conventional normalization with [co, cj] = 1 and CqCo 
determining the operator of the photon number. 

A more realistic field consists of a certain number of 
modes propagating within a paraxial beam of a cross- 
sectional area S±. If we keep the same notation cq for 
the field operators describing the field amplitude in the 
beam, their commutator becomes [co,Cq] = Aj/V where 
Aj is the number of modes. The total photon flux density 
in the beam is then given by Q = v gr S±(^F\cQCo\^F)- 
It is convenient to go from a discrete set of modes to a 
continuous spectral interval Aw <C u>. The density of 
states in a volume V is equal to -q — Vk 2 /8n 3 v gr and the 
wave vectors of the modes constituting a beam occupy 
the solid angle Ao « 4Tr 2 /k 2 S±. As a result, we arrive 
at the following commutation relations for the operator 
of the field amplitude and its spectral harmonics: 
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(4) 
(5) 
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Here the spectral decomposition of the field amplitude 
operator is defined as Co = J A c.Q U e~ lu)t dui . 

The equation of motion for the field amplitude opera- 
tor of each of the two signal fields has the same form as 
the wave equation for a classical field amplitude: 



d_ 
dt 



d 



E d(0J 2 n 2 )/duj 
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(7) 



Equation ([7| includes all the relevant effects: linear 
dispersion determines the group velocity of the wave, 



whereas the slowly varying polarization amplitude Po on 
the right-hand side includes nonlinearity, dissipation, and 
fluctuations. At the boundary z b between the medium 
and the vacuum, the boundary condition for the field 
operator takes the form 
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which satisfies the conservation of the Poynting flux. 
Eqs. ([3]) are to be solved together with Eq. Q for both 
signal fields in order to determine the generated signal 
and noise. 

In the four- wave mixing process depicted in Fig. 2, the 
total field consists of the four waves, 

E = e {+) E HF e~ lUHFt+lkHFZ + e {+) E LF e- lujLFt+lkLFZ 



+ c.c. + e (+) ^oC(+)e~ 2 "<+» t+lfc( + )Z 
+ e ( _)£:oC(_)e~ l "(-> t+4fe <-> 2 + h.c. 



(9) 



in which two strong classical pump fields at high and 
low frequencies (denoted as HF and LF) are resonant to 
the corresponding transitions between the Landau levels, 
ujhf = W41 and lulf = ^21, whereas two signal fields are 
described by operators and their frequencies may have 
a detuning, W( +! _) = 0^43^2 =F A, A <C satisfying 
the frequency-matching condition uihf = ^lf + L0 {+) + 
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We also assumed that ui 
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(uj) in the 



normalization constant Eq. 

Solving the density-matrix equations ^ in the steady 
state and in linear approximation with respect to weak 
signal fields is reduced to a straightforward algebra. Op- 
timal conditions for the entanglement are realized when 
the Fermi level is close to the state 1 (n = -2) and the 
populations of all states above are low. This is possible 
when the magnetic field is strong enough, fc^T < hui c , 
and Rabi frequencies of the pump fields are below satu- 
ration: \Qhf,lf I (7)- Here the Rabi frequencies are 



defined as £Ihf 



d^Ejjp 



assume for simplicity that all scattering rates j mn are 
of the same value (7). The latter assumption can be 
easily dropped once the relaxation rates are known for 
any particular sample. If, in addition, the detuning is 
sufficiently large, (7) <§; A, the populations of the ex- 
cited states are mostly due to the Langevin noise terms 
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and we 



F(+ t _) = -£32,43 in Eqs. (|3j). Solving the density-matrix 
equations in the steady state and neglecting the terms of 
the order of ((7) /A) 2 , we arrive at the following expres- 
sion for the operator of the polarization amplitude at the 
frequency of the signal fields: 
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and we denoted fi 2 = Hhf^lf an( l d = HevF/w32- 

Using the polarization ( 10 )) as a source in Eqs. 0, we 
obtain the following coupled equations for the signal field 
operators: 
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where the coefficient of the parametric coupling is 

/,.,\2 



k = 2wx 



c 2 (k) 



and the noise term 



G (+--) 



(w) 2 rf (+ _)F (+ _ } 
c 2 {k) E A 



(13) 



Here we again neglected a small difference between the 
central frequencies of the signal fields in the pre-factors, 
assuming (*)t+) = = (u>) and (fc) = (lo)/c. 

In the optimal limit of \Cl p \ <C (7) <C |A|, the noise 
terms and the Raman scattering of the pump fields into 
the signal modes can be neglected and the solution for the 
fields exiting a layer of thickness L takes a particularly 
simple and transparent form: 

c (+) (L,t) = cosh(r)c (+) (0,t - L/v gr ) 

- ie ie sinh^cK (0, t - L/v gr ) 
£(_)(£, t) = cosh(r)c(_) (0,t-L/v gr ) 

- ie ie sinh(r)cJ +) (0, t - L/v gr ) . (14) 

Here the parametric gain factor r = \k\L and k = \ii\e . 
Eqs. (14) clearly show the emergence of quantum cor- 



relations between the signal photons with opposite cir- 
cular polarizations. In particular, consider the boundary 
condition at z — corresponding to a completely uncor- 
rected state of vacuum fluctuations within the spectral 
bandwidth Au. Then one can obtain from Eq. ( 14 1 that 



the photon fluxes in two signal fields exiting the layer at 
z = L are completely correlated: 



(0\Q {+) (L)\0) = (|Q ( _)(i)|0) 

(Q|((5 (+) (L)-<3 ( _ ) (L)j |0)=0 



sinh 2 r, 

2 2lT (15) 



Here Q( +) _)(L) = cSxcj + _JL)c( + _^(L) are operators 
of the photon fluxes. The correlated (+) and (-) photons 
can then be used to prepare the desired polarization- 
entangled states, e.g. the analog of the Bell states pQ. 
The second equation in ( 15 ) corresponds to the Manley- 



Rowe relations for the parametric process. It also follows 
from Eq. ( 14 ) that the scheme could be used to amplify 
the light with a nonclassical statistics or exchange the 
statistical properties between (+) and (-) photons. The 



magnitude of Aw is likely to be limited by the bandwidth 
of a detection system. 



If noise terms G^ in Eq. ( 12 ) are taken into account, 



the field equations are still straightforward to solve, al- 
though the solution becomes more cumbersome. As a 
result, the photon fluxes in Eq. (15) acquire additional 
noise terms: 
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are of the order of the parametric coupling terms: 
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From this solution one can see that the noise contribu- 
tion can be neglected if |A| 3> (7) provided the paramet- 
ric gain is high enough: r > 1. For a weak amplification 
t <C 1 the condition for a large signal to noise ratio is 
more stringent: A ^> (7)/t. If this condition is not 
satisfied or if one of the states 2, 3, or 4 acquires a large 
population, the noise is always comparable to or greater 
than the signal in the steady state. Then the entangled 
photons can be generated only in the pulsed regime dur- 
ing the time of the order of a few relaxation times I/7. 
This is usually the case in resonant schemes of entangled 
photon generation in atomic vapors [TH1 [20] ■ 
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FIG. 3: Parametric gain r per monolayer of graphene as a 
function of normalized detuning of the signal fields A / (7) for 
the pump field intensity |fi p | 2 = 0.1(7) 2 . 
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FIG. 4: Parametric gain r per monolayer of graphene as a 
function of normalized strength of the pump field |fi p |/(7) for 
the detuning A = 10(7). 

The above analytic results were derived in the limit of 
|fi p | <C (7) <C A. In the general case the equations can 
be solved numerically, including the effects of the optical 
pumping of electrons to excited states and optical satura- 
tion. The resulting parametric gain r per one monolayer 
of graphene is plotted in Figs. 3 and 4 as a function of 
the frequency detuning and the pump intensity. As seen 
from the figures, the magnitude of r is around 0.01 for 
A ~ IO7 ~ 100Op. This corresponds to a photon flux of 
about 10~ 4 Aw/27r. To increase the value of r for a higher 
rate of the twin photon generation, one can use a stack of 
graphene monolayers or a thin layer of graphite. Recent 
studies have demonstrated that a graphite layer consist- 
ing of ~ 100 monolayers maintains high carrier mobility 
and graphene-like Landau levels near the H-point of the 



Brillouin zone of graphite [131 Q3] . The optimal thick- 
ness is a tradeoff between the pump absorption and the 
desired output photon flux. 

We showed that graphene placed in a magnetic field 
can serve as an efficient, tunable source of polarization- 
entangled photons in the mid-infrared and THz frequency 
range. The proposed scheme can operate at high temper- 
atures if the magnetic field is high enough to prevent ther- 
mal excitations, fc^T <C huj c . A thin layer of graphene 
can be easily integrated with semiconductor laser chips 
and optoelectronic circuits to make a compact setup. 
Similar mechanism of entangled photon generation could 
exist in topological insulators where the surface states 
have a massless dispersion and demonstrate a similar pat- 
tern of Landau levels oc V Bn and magneto-absorption 
|21) . Interestingly, the band velocity vp ~ 8.5 x 10 7 cm/s 
for surface states in Bio.91Sbo.09 inferred from measure- 
ments in |21j is very close to the one in graphene, which 
suggests an optical nonlinearity of similar strength. Wc 
hope that our results will stimulate active experimental 
research in nonlinear and quantum optics of graphene 
and other materials with similar electronic properties. 
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